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Preface 


The effectiveness of our nuclear arsenal to act as a nuclear 
deterrent hinges in part upon the credibility of its 
survivability. Simply stated, if a potential adversary does 
not believe our arsenal would survive deployment during a 
nuclear exchange then our weapons present a reduced threat 
to him. To ensure the survivability of our nuclear 
weaponry, we must be able to test critical systems and 
components in environments that realistically simulate the 
stresses created by a nuclear environment. Currently, there 
are serious gaps in our weapons effects simulations 
capabilities; should a comprehensive nuclear test ban treaty 
be signed, these gaps would become even more severe. 

This work examined the feasibility of using the hard X-ray 
pulse generated in the proposed Laboratory Microfusion 
Facility (LMF) for hard X-ray weapons effects simulation. 

The construction of the LMF is a necessary next step in the 
development of inertial confinement fusion power for 
commercial electricity generation; the application of the 
LMF to weapons effects testing could significantly enhance 
the national security. The increased funding resulting from 
weapons effects applications would encourage the development 
and construction of the LMF and thereby encourage the 
development of nuclear power for peaceful purposes. 



There are a number of people who provided help and valuable 
guidance to me on this project. Captain Hike Tobin, my 
sponsor from LLNL, provided much of the background 
information used in this work and his enthusiasm for the 
project was contagious. Major Denis Beller, my advisor from 
AFIT, was always willing to talk with me when I arrived 
(usually without an appointment) at his door, bursting to 
share something new I had learned. He was also especially 
helpful in sorting out problems with the computer code and 
in providing general guidance and direction to the work. My 
wife, Kathy, and my children, Jason, Mark and Hannah, 
supported me through what was a difficult time for all of 
us; I cherish them deeply and am thankful for such a 
wonderful family. More than anyone else, however, I would 
like to express my gratitude to my Lord Jesus Christ, who is 
my breath, and who was the real reason I came to AFIT in the 
first place. 


Don F. Nichols 
March 6, 1989 
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Abstract 


This work examined the feasibility of using the proposed 
Laboratory Microfusion Facility (LMF) to simulate the 
effects of hard X rays from a nuclear weapon. The LMF will 
be an inertial confinement fusion facility for the testing 
of high-gain deuterium-tritium (DT) pellets, and will 
produce a pulse of hard X rays and neutrons over a very 
short time interval. The purpose of the study was to 
determine (1) the intensity that remains after the X-ray 
pulse is stretched to the 10's of ns required for weapons 
effects work, (2) the effects of the temporal stretching on 
the final X-ray spectrum, and (3) the intensity of the 
neutron fluence and neutron-generated gamma rays. Using the 
MORSE-CG Monte Carlo code, a spherical shell of 6 LiH with a 
thickness of 60 cm was found to spread an instantaneous 10 
keV X-ray pulse temporally over a period of 16 ns. The 
spectrum of the resultant X-ray pulse was somewhat softer 
than that of the input spectrum; the 10 keV input X-ray 
spectrum was peaked at 37.5 keV while the final spectrum was 
peaked at 25 keV. Eighty-seven percent of the X-ray energy 
fluence was lost to the LiH, resulting in a final fluence of 
10.8 cal/cm 2 and a final dose of 9.3 Mrads(Si) 250 cm from 
the source. Fifty percent of the X-ray dose arrived during 
an interval of 6 ns, giving a dose rate of 8*10 14 
rads(Si)/sec. The only serious problem with the use of the 
LMF for hard X-ray simulation was the fluence of one-MeV- 
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equivalent neutrons. The minimum fluence was on the order 
of 1*10 14 n/cm 2 , but preliminary work has indicated that 
this value could be reduced by an appropriate shadow shield. 
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SPECTRAL AND TEMPORAL FIDELITY OF A 
HARD X-RAY WEAPONS EFFECT8 SIMULATION TEST 
IN A HIGH-GAIN ICP FACILITY 


I. Introduction 


Background 

In the design of weapons systems it is essential that 
the weapons designer be able to assess not only the 
lethality of his weapon, but also its vulnerability to 
real-world threat environments. A weapon is only as 
good as its ability to survive its journey to the 
target; if, for example, a re-entry vehicle is rendered 
impotent by the effects of a precursor nuclear burst, it 
has been launched in vain. For this reason, an entire 
scientific community has evolved around the work of 
weapons vulnerability, lethality and effects (VLE). 

This community has spent considerable time and effort 
developing simulators that attempt to re-create the 
radiation effects which weapons systems would experience 
during a nuclear exchange. Such simulation is quite a 
task, as the nuclear environment produces a complex set 
of synergistic phenomena through the interaction of X 
rays, gamma rays, neutrons and weapons debris with the 
weapon system. Adequate simulation of the nuclear 









environment requires special attention to the spectral 
and temporal characteristics of the actual stockpile- 
to-target threat sequence (STS); the degree to which a 
particular simulator duplicates the stresses produced by 
these characteristics is known as the fidelity of the 
simulator. 

While an examination of all the effects of a nuclear 
detonation is outside the scope of this paper, the 
problem of weapons effects simulation can be illuminated 
by a brief examination of some of the principal photon 
and neutron effects that must be considered in weapons 
effects simulation. 

Photon Effects , The X-ray flux from an endo- 
atmospheric nuclear burst is quickly absorbed by the 
atmosphere, and its effect on other weapons systems is 
usually only of minor importance. The X-ray flux from 
an exo-atmospheric nuclear burst is diminished only by 
spherical divergence, however, and large X-ray fluences 
can be present at great distances from the burst. These 
X rays can seriously affect vehicles, such as satellites 
and re-entry vehicles (RV's), that move outside the 
earth's atmosphere. High fluences of low-energy X rays 
on the surface of a RV can ablate its heat shield. If 
the shield is sufficiently damaged it may not provide 
adequate protection and the RV may burn up during re- 
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entry into the atmosphere. Large doses of low-energy X 
rays over short time scales can give rise to a thermo- 
mechanical-shock response that propagates through the 
structural material and destroys the vehicle by 
spallation of internal surfaces(1:5.1,22). X and gamma 
rays can also dislodge electrons by the photoelectric 
effect. These electrons, if generated over a 
sufficiently short time and in sufficient quantity, can 
give rise to electric currents that interact to create a 
system generated electromagnetic pulse (SGEMP). This 
pulse can disrupt the internal electronics and lead to 
total system failure. Similarly, hard X and gamma rays 
can interact directly with the electronics and generate 
electrons through photoelectric and Compton 
interactions, and by pair production. In the 
semiconductors used in modern electronics, the electrons 
knocked free and the "holes" where they once resided can 
seriously change the characteristics of the affected 
device, causing changes in device states (upset) and 
voltage thresholds. High dose rates can result in the 
generation of large photocurrents that can destroy a 
device or burn out an entire circuit(2:275-296). 

Neutron Effects . Neutrons can interact with the 
fissile material in a nuclear warhead and generate 
enough heat to asymmetrically soften its components, 
preventing a symmetric final compression and causing the 






warhead to fail. Neutrons can also interact with the 
lattice structure of the semiconductors used in modern 
electronics to create defects that increase the 
resistivity of the semiconductor and reduce the minority 
carrier lifetime. This changes the device operating 
characteristics and may render affected devices and the 
components of which they are a part inoperable(3:385). 

The neutron and photon effects presented above generate 
complex interactions between electronic devices at both 
the component and the system level. For this reason, it 
is desirable to test the effects of a species of 
radiation on an entire system, not just its components. 
Current above-ground test (AGT) X-ray simulation 
facilities are not adequate to test large items such as 
re-entry vehicles and satellites. Further, 
simultaneously achieving high fidelity with respect to 
X-ray spectrum, dose, dose rate and the rate with which 
the dose rate changes is quite difficult in AGT 1 s and, 
perhaps surprisingly, even in underground tests (UGT's) 
using actual nuclear devices. 

In view of this difficulty, the development of high- 
fidelity X-ray weapons effects simulators is of great 
importance to assuring the survivability of our nuclear 
weaponry and, therefore, of great importance to our 
national security. One facility currently being 









evaluated for its weapons effects simulation 
capabilities is the Laboratory Microfusion Facility 
(LMF), a proposed facility to be used to test and 
develop high-gain inertial confinement fusion (ICF) 
pellet designs. Recent advances in the ICF field have 
led to the possibility of ICF pellets that could 
generate energy gains of 100 fold. The most promising 
pellet designs are those that use deuterium and tritium 
(DT) as the fusion fuel. These high-gain designs are 
predicted to be able to produce an energy pulse of 1 GJ 
for a 10 MJ driver at a power of 10 17 to 10 18 watts. 

This pulse will consist of a large neutron component and 
an X-ray component that, if it can be suitably tailored, 
could be used to simulate the hard X-ray pulse (X-ray 
energies greater than 10 keV) a device or weapons system 
might experience in space from a nuclear burst. The 
size of the facility and the amount of hard X-ray energy 
released by the fusion of the ICF pellet could make the 
LMF suitable for the hard X-ray testing of entire 
systems such as satellites or re-entry vehicles(4:5). 
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Problem and Scope 


Several issues arise in regard to tailoring an ICF X-ray 
pulse for use in hard X-ray VLE simulation. The first 
issue to be resolved is the degree to which a pulse of X 
rays can be stretched temporally without sacrificing an 
acceptable dose outside the scattering medium. The ICF 
pulse is very short but Bridgman suggests that an exo- 
atmospheric X-ray pulse would be as long as several tens 
of nanoseconds(1:5.22). Substantial temporal stretching 
of the initial X-ray pulse would be required to achieve 
dose-rate fidelity; in the process the pulse could be 
substantially attenuated. To be attractive for hard X- 
ray simulation a dose greater than 1 Mrad(Si) must 
remain after the pulse has been lengthened. The second 
issue to be considered is the effect of the scattering 
process on the X-ray spectrum. In most materials, the 
absorption cross section for low-energy photons is much 
greater than for high-energy photons; the spectrum of an 
X-ray pulse passing through a material is typically 
hardened. Excessive pulse hardening could render the 
LMF unsuitable for VLE work. The final issue to be 
addressed is the neutron and neutron-generated gamma-ray 
fluence outside the scattering shell. High-energy 
neutrons can be a source of radiation damage to both 
structures and electronics. Neutron fluences greater 
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than 10 13 MeV-equivalent neutrons/cm 2 would interfere 
with the X-ray effects being studied. The effects of a 
pulse of high-energy neutron-generated gamma rays could 
overwhelm the X-ray effects and deliver an unacceptably 
high total dose to electronics. 

To resolve the issues raised above, the MORSE-CG Monte 
Carlo computer code was used to calculate hard X-ray and 
neutron doses outside a lithium hydride (LiH) shell in 
an LMF chamber. Lithium hydride was chosen as a 
scattering medium because of it's high scatter-to- 
absorption cross-section ratio for photons. The problem 
approach is presented in the following section. 

Approach 

The effects of various blankets of LiH on the ICF pulse 
were studied using the Multi-group Oak Ridge Stochastic 
Experiment - Combinatorial Geometry (MORSE-CG) computer 
code, distributed by the Radiation Shielding Information 
Center (RSIC) at Oak Ridge National Laboratory (ORNL). 
The steps to the investigation were as follows: 


1. Compile and run the sample problems and code as 
supplied by RSIC, and fix any bugs discovered. 







2. Premix cross-section data required by the 
specific problem, using the MORSE module 


XCHECKER. 

3. Develop a simple geometry (using the MORSE module 
PICTURE) with which to check user-written 
routines. 

4. Develop user-written routines to provide the 
specific source and detector responses required 
by the specific problem of interest. 


5. Test the user-written routines using single group 
input spectrum and simple geometries. 



6 . 


Develop the specific geometry to be used for the 
actual problem and test it thoroughly using the 


MORSE module PICTURE. 


7. Develop input spectral data for both X rays and 
neutrons, gamma-generation information (from the 
XCHECKER run), and any KERMA (Kinetic Energy 
Released in Material) factors desired. Combine 
these data with the geometry data to create the 
MORSE input file. 


8 . 



Run the problem and examine the results. Alter 
the geometry and rerun the problem to investigate 
relevant phenomena. 
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Of the above steps, steps one, three, four and five were 
actually preparatory steps, and are not the subject of 
this thesis. They did, however, consume an enormous 
amount of time. The routines developed in step four and 
the data developed in the remaining steps will be 
discussed in full in the body and appendices of this 
report. 

Organization of Presentation 

A detailed description of this investigation is 
presented in the chapters that follow. Chapter 2 
presents background information on the MORSE code and 
briefly discusses the theory on which it operates. 
Chapter 2 also presents a brief discussion of the user- 
written routines and the geometry, cross-section and 
KERMA data. A large number of problems were run and 
analyzed in this work; Chapter 3 presents several 
selected problems, the results and discussion of which 
are illustrative of the major findings of this work. 

For clarity, each problem is discussed as it is 
presented. Chapter 4 presents the conclusions drawn 
from this work and gives specific recommendations for 
future work. A discussion of the functions of ^ach of 
the user-written routines is presented along with their 








source code in Appendix A end, finally, a sample data 
file is presented in Appendix B. 
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II. The MORSE Code 


Description off the Code 

Theory and Operation . The Multi-group Oak Ridge 
Stochastic Experiment - Combinatorial Geometry (MORSE- 
CG) computer code is a statistical approach to solving 
the Boltzman transport equation for both photons and 
neutrons. The code is based on an interpretation of the 
macroscopic cross section in which the cross section for 
a given nuclear reaction represents the probability per 
unit path length of that reaction occurring. A thorough 
discussion of the mechanics of the code would consume 
many volumes and is beyond the scope of this paper; a 
brief overview of the code is presented, however, in the 
following paragraphs. 

The MORSE code is a multi-group, Monte Carlo code; the 
particle energy distribution is divided into user- 
defined groups (or bins), and the angular distribution 
is represented by an expansion of Legendre polynomials 
up to the sixteenth order (only a third order (P3) 
expansion was used in the calculations performed for 
this work). A simple text file is used to enter 
specific problem data. This file contains the multi¬ 
group source spectrum, its energy structure, the problem 
geometry and cross-section information. The source 
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spectrum is simply a listing of the number of source 
particles in each of the user-defined bins. For a 
coupled, neutron-photon calculation, the upper bins 
(which begin with bin number 1) of the source spectrum 
contain the neutron data and the lower bins contain the 
photon data. The code reads the spectral data from the 
input file, normalizes the total number of particles in 
the spectrum to 1, and uses the normalized spectrum to 
generate a batch of particles with the same spectral 
composition as the input spectrum. Each of these 
particles is assigned an initial direction, weight and 
velocity, based on the energy of the particle and 
various other parameters supplied by the user. Each 
particle is then "tracked" as it travels through and 
collides with the media which comprise the problem 
geometry. 

At each collision site, particles can undergo a number 
of events. Particles can scatter into a new direction 
and/or into a new energy group, and may generate new 
particles through such mechanisms as pair production, 
neutron capture, or (n,2n) reactions. MORSE determines 
what happens at each collision through the use of cross 
sections, probabilities and random numbers. MORSE 
determines the chronological age of each particle by 
maintaining a running sum of the distance between 
collisions and boundary crossings, divided by the 
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particle speed. Particles are never absorbed in MORSE; 
at each collision site the pre-collision weight of the 
particle is reduced by the ratio of the cross section 
for absorption to the total cross section. This 
prevents absorption from decreasing the number of 
particles in the system and reduces the average 
contribution that single particles make to the total 
fluence at the detector, hence improving the statistics. 
Particles are tracked through the system until their 
weights become very small or they escape across the 
outer boundary of the geometry. 

Eventually, a particle's weight becomes so small that 
the particle can no longer make a significant 
contribution to the total fluence. To continue to track 
such a particle only uses up computer time and does not 
move any closer to the solution of the original problem. 
To avoid tracking numerous particles with tiny weights, 
methods are available for 'combining' several low-weight 
particles to form a single particle, the weight of which 
equals the combined weights of the particles eliminated. 
This work used a method known as "Russian roulette". In 
Russian roulette, when a collided particle's weight 
falls below a user specified weight cutoff, the ratio of 
the current particle weight to a user-supplied survival 
weight is compared to a random number generated by the 
computer. If the random number is greater than the 
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ratio, the particle is killed and its history is 
terminated. If the random number is less than or equal 
to the ratio, the particle is given the user-supplied 
survival weight and allowed to 
continue(5:4-1,2,4-10,11). 

The User-written Routines 

The MORSE code is a system of subroutines called by a 
main program module. Standard subroutines included in 
the package perform functions that are common to all 
applications (such as the random walk process and 
statistical bookkeeping), but the user is expected to 
generate problem-specific routines. The routines to be 
generated for this work were the main routine (that 
calls MORSE), and subroutines named Bankr, Direc, Gtmed, 
Relcol, Sgam, Sdata, and Source. These routines were 
adapted from routines that already existed; the 
functions of these routines and their source code are 
presented in Appendix A. 
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Input Geometry Descriptiona 

An aluminum shell, 1 meter in thickness and with an 
inside radius of 5 meters, was used to simulate the LMF 
chamber. This shell, with it s contents, is illustrated 
in Figure 2.1, which is a compressed output from the 
MOPCZ module PICTURE. The chamber is surrounded by free 
space. 

A pyrolytic-graphite cone with a height of two meters 
and a base radius of twenty-five centimeters was used to 
simulate a generic RV. It is illustrated in Figure 2.2. 
It is mounted to the wall of the chamber by two aluminum 
support cylinders, each five centimeters in radius. To 
maximize the uniformity of the X-ray fluence on the edge 
of the graphite cone nearest the source, the cone is 
positioned with the surface closest to the center of the 
chamber parallel to the vertical ( or 'z') axis of the 
chamber (see Figure 2.1) and such that the 'x' axis of 
the chamber intersects the midpoint of the cone. To 
simulate semiconductor material in an aluminum casing 
inside an RV, a spherical silicon shell with an outer 
radius of 2.5 cm, an inner radius of 1.1 cm and a 0.254 
cm casing of aluminum is located 1.25 cm inside the 
surface of the graphite cone closest to the source. 
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Note: Letters represent materials; g is graphite, L is 
Lithium hydride, W and # are aluminum, and . and 0 are 
free space. The chamber radius is 5 meters. 

Figure 2.1: Basic Target Chamber Geometry 
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Note: Letters represent materials; g is graphite, # is 
aluminum, I is ice and . represents free space. The RV 
is 2 m in height with a base diameter of 50 cm. 

Figure 2.2: Generic Conical Re-entry Vehicle 


A shell of lithium hydride (LiH) was chosen to 
temporally stretch the X-ray pulse length and attenuate 
the neutron fluence. The PICTURE depiction of this 
shell is presented in Figure 2.3. The shell is located 
at the center of the LMF chamber, concentric to the 


























































































































































location of the source. The inner radius of the basic 
shell is 40 cm and the outer radius is 80 cm. It is 
supported by an aluminum pedestal five centimeters in 
radius. The ICF pellet would initially be supported by 
a stem; this was modeled using a tapered aluminum stem 
which extends through the LiH shell from the pedestal 
and ends five centimeters from the center of the 
chamber. 
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Note: Letters represent materials; L is lithium 
hydride, § is aluminum, i is ice and . represents 
free space. The inner radius is 40 cm and the 
outer radius is 80 cm. 


Figure 2.3: The LiH Moderating Blanket 


All surfaces outside the LiH shell and inside the 
aluminum chamber are coated with a layer of low density 
ice with a thickness of 1 cm. Each structure, such as 
the graphite cone, the wall, or the LiH shell can be 
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'turned off' for a given problem by setting its medium 
number in the data input file to 1000, the number 
signifying an internal void. A structure is 'turned on' 
by setting its medium number to the number representing 
the cross sections to be used for the structure. In 
this fashion the effects of the various chamber 
components can be examined separately or in combination. 

The Input Spectrum 

Since the actual output spectrum from an ICF capsule is 
classified, an actual spectrum was not used in this 
work. However, Chen(6:14) states that for a deuterium- 
tritium (DT) reaction, a plasma is required with a 
temperature in the 10 keV range. Thus, the X-ray 
spectrum used in this work is a 10 keV Planckian 
blackbody distribution. The distribution of neutrons 
with energy is presented in Table 2.1 and shown 
graphically in Figure 2.4. According to Friedlander 
et. al.(7:546), the neutrons produced in the DT reaction 
carry an average of 14.1 MeV of the 17.6 MeV released in 
the reaction, so the neutron spectrum used in this work 
is composed of predominantly 14 MeV neutrons with a 
lower-energy tail to account for neutrons which collide 
in the DT pellet. 
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Energy Range 

Neutrons in Range 

(MeV) 

(% Total) 

14-10 

49.8 

10-1 

39.3 

<1 

10.9 

















Cross-section Data 


The microscopic cross sections used in this work were 
supplied by RSIC at ORNL as a part of the FEWG1 RSIC 
Data Library Collection^: 1) . This library supplies 
cross-section data for neutrons with energies from 10 
micro-eV to 19.6 MeV and for photons with energies from 
10 keV to 14 MeV. This range of photon cross sections 
limited this work to calculations with hard X rays. The 
material densities used to generate the macroscopic 
cross-sections are listed in Table 2.2. 


Table 2.2: Materials and Densities 


Material 

Density 

gm/cc 

atoms/b-cm 

Aluminum 

2.7 

6.022E-2 

LiH 

0.6 

4.520E-2 

Graphite 

2.2 

1.104E-2 

Ice 

0.1 


Hydrogen 


3.346E-3 

Oxygen 


6.691E-3 

Boron 

2.30 

1.281E-1 

In LiH 

0.23 

1.281E-2 

Silicon 

2.33 

4.996E-2 

Lead 

11.3 

3.284E-2 
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Detector Response Functions 

The basic MORSE code produces fluence output in terms of 
particles per square centimeter per source particle. 
Because energy fluerces were desired, the basic code was 
modified to give the fluence in terms of 
eV/cm 2 /source particle. In many situations, however, it 
is not the raw fluence that is of interest, but rather 
the energy deposited in the material. Conveniently, the 
FEWG1 cross-section library used in this work also 
includes a set of KERMA (Kinetic Energy Released in 
MAterial) factors. The units on these factors are 
(rads-barn-gram)/(particle-atom) where 'particle' refers 
to the incident particle (neutron or photon) and 'atom' 
refers to the media on which the particle is impingent. 
To convert these units to rads/(eV/cm 2 ) it is necessary 
to divide the KERMA factor for each energy bin by the 
average bin energy in eV and by the atomic weight in 
grams and to then multiply by 10 -24 cm 2 /barn. 

The factors discussed above are convenient for X-ray and 
gamma-ray doses, but neutron fluences are typically 
reported in terms of one-MeV-equivalent neutrons/cm 2 
(2:379,380). For a given fluence of neutrons with 
energy in dE about E, the one-MeV-equivalent fluence is 
defied as the fluence multiplied by the ratio of the 
KERMA factor for that energy divided by the KERMA factor 
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for a 1 MeV neutron, which Messenger gives for silicon 
as 1.52*10 -9 Rads-barn-gram/particle-atom. 


Using these values, response functions were calculated 
for use with the MORSE code. These function values are 
presented in Tables 2.3 and 2.4. 

Table 2.3: Photon Detector Response Functions 



Upper 

Energy 

Boundary 

(eV) 


Output Units per Source Particle 


eV/cnU 


rads/C) 

(* 10 * 7 ) 


rads(Si) 
(* 10 17 ) 


Energy 


Response Function Values 


+ 1.4000 

+7 

1.0000 

+7 

8.0000 

+ 6 

7.0000 

+ 6 

6.0000 

+ 6 

5.0000 

+ 6 

4.0000 

+ 6 

3.0000 

+ 6 

2.5000 

+ 6 

2.0000 

+ 6 

1.5000 

+ 6 

1.0000 

+ 6 

7.0000 

+ 5 

4.5000 

+ 5 

3.0000 

+ 5 

1.5000 

+ 5 

1.0000 

+ 5 

7.0000 

+ 4 

4.5000 

+ 4 

3.0000 

+ 4 

2.0000 

+ 4 



2.190 + 

30.12 

2.331 

29.76 

2.435 

29.64 

2.540 

29.91 

2.682 

30.56 

2.870 

31.52 

3.126 

33.01 

3.392 

34.85 

3.624 

36.69 

3.914 

39.22 

4.267 

42.59 

4.581 

45.74 

4.733 

47.40 

4.698 

47.53 

4.369 

46.92 

3.698 

55.79 

3.390 

92.25 

3.453 

241.9 

10.43 

924.2 

16.35 

3305.0 

92.77 

17630.0 
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Table 2.4: Neutron Detector Response Functions 


Upper 

Energy 

Boundary 

Output Units per Source Particle 

* / 2 
n /ever 

rads(C) 

rads(Si) 

Energy(eV) 

Response Function Values 

1.9600 

+ 7 + 

3.208 -06 

1.216 -17 

1.047 -16 

1.6900 

+7 

3.429 -06 

1.186 -17 

1.119 -16 

1.4900 

+7 

3.543 -06 

1.131 -17 

1.156 -16 

1.4200 

+7 

3.553 -06 

1.167 -17 

1.16 -16 

1.3800 

+7 

3.490 -06 

1.156 -17 

1.139 -16 

1.2800 

+7 

3.396 -06 

1.182 -17 

1.108 -16 

1.2200 

+7 

3.254 -06 

1.100 -17 

1.062 -16 

1.1100 

+7 

3.059 -06 

1.130 -17 

9.983 -17 

1.0000 

+7 

2.828 -06 

1.549 -17 

9.229 -17 

9.0500 

4-6 

2.831 -06 

9.79 -18 

9.239 -17 

8.1900 

+ 6 

3.070 -06 

1.506 -17 

1.002 -16 

7.4100 

+ 6 

2.093 -06 

6.93 -18 

6.831 -17 

6.3800 

+ 6 

1.155 -06 

1.055 -17 

3.770 -17 

4.9700 

+6 

1.156 -06 

1.068 -17 

3.771 -17 

4.7200 

+ 6 

9.899 -07 

1.442 -17 

3.230 -17 

4.0700 

+ 6 

9.161 -07 

2.472 -17 

2.990 -17 

3.0100 

+ 6 

1.136 -06 

2.212 -17 

3.706 -17 

2.3900 

+ 6 

1.082 -06 

1.830 -17 

3.532 -17 

2.3100 

+ 6 

1.403 -06 

1.954 -17 

4.579 -17 

1.8300 

+ 6 

1.362 -06 

2.180 -17 

4.435 -17 

1.1100 

+ 6 

1.838 -06 

2.903 -17 

5.999 -17 

5.5000 

+ 5 

2.723 -06 

3.643 -17 

8.903 -17 

1.5800 

+5 

4.078 -07 

4.794 -17 

1.331 -17 

1.1100 

+ 5 

1.051 -06 

4.844 -17 

3.428 -17 

5.2500 

+ 4 

0.000 

5.019 -17 

3.190 -17 

2.4800 

+4 

0.000 

5.313 -17 

3.387 -17 

2.1900 

+4 

0.000 

5.121 -17 

3.661 -17 

1.0300 

+4 

0.000 

4.906 -17 

4.367 -17 

3.3500 

+ 3 

0.000 

5.018 -17 

4.585 -17 

1.2300 

+ 3 

0.000 

5.196 -17 

4.778 -17 

5.8300 

+ 2 

0.000 

4.370 -17 

4.217 -17 

1.0100 

+ 2 

0.000 

4.841 -17 

7.061 -17 

2.9000 

+ 1 

0.000 

4.674 -17 

1.984 -16 

1.0700 

+ 1 

0.000 

5.090 -17 

8.182 -16 

3.0600 

+ 0 

0.000 

6.606 -17 

4.488 -15 

1.1300 

+ 0 

0.000 

1.208 -16 

1.994 -14 

4.1400 

-1 

0.000 

1.093 -15 

3.054 -13 


+ Read as 1.96*10 7 eV * one-Mev-equivalent 
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III. Results 


Organization 

This chapter presents a selection of the results of this 
investigation. These results illustrate the significant 
information gained from the rather large number of 
problems run. For clarity, the problems presented are 
arranged by the nature of the information they contain 
and not necessarily in the order in which they were run. 
Further, to provide a clear and logical development of 
the conclusions of this work, each problem is discussed 
as the data pertaining to that problem are presented. 

Dolan suggests an X ray energy fraction of 2 to 15 
percent of a total yield of 30 to 3000 MJ(9:515). This 
work used a total energy output (ignoring debris) of 
1000 MJ, and allocated 867 MJ as neutron source energy 
and 133 MJ as X-ray energy. This arbitrary energy 
distribution is used to give the order of magnitude of 
the effects being discussed. 

The results are presented according to outline presented 
on the following page. 
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X-ray Results (source neutrons off) 

LiH Shell in an LMF Chamber 

Neutron Results (source X rays off) 

LiH Shell in an LMF Chamber 

Combined (source X rays plus neutrons) Results 

LiH Shell in an LMF Chamber 

Dose in Rads(Carbon) and Rads(Silicon) 
at a Graphite Cone 

Reduction of Neutron and Gamma-ray Fluence 
Addition of Boron-10 to LiH 
Shielding of Wall Pulse 

Substitution of 6 LiH and Realistic Walls 
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X-ray Fluence, Dose and Arrival Tinea 


The first issue of concern is the effect of the LiH 
scatterer on the total fluence and dose of hard X rays. 
To resolve this issue, a series of calculations was 
performed using only the X-ray portion of the input 
spectrum (i.e. with the neutrons turned 'off'). The 
temporal effects were restricted to those of the LiH 
shell and the chamber wall by turning all other material 
in the problem "off". Since the time interval over 
which the actual ICF X-ray pulse would be emitted from 
the DT capsule is very short compared to the length to 
which it must be stretched, the X-ray source was modeled 
as an instantaneous point source located at the center 
of the chamber. The X-ray energy fluence and dose as a 
function of time and LiH shell thickness are presented 
graphically in Figure 3.1 and numerically in Table 3.1. 
Note that although the information presented in Figure 
3.1 and in the figures to follow is presented in the 
form of line graphs, the data are more accurately 
displayed in the form of a histogram. Each point 
actually represents the midpoint value of an energy or 
time bin; line graphs are used to clarify the 
interpretation of overlaid plots. 
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Table 3.1: Variation in X-R^y Fluence, Dose and Arrival 
Interval with 'LiH Thickness 


Outer 

radius 

(cm) 

LiH 

Shell 

(cm) 

Fluence 

(cal/cm 2 ) 

Dose 

MRads(Si) 

T50 

(ns) 

T9 0 
(ns) 

40 

PI 

80 (1) 

72 (2) 

<1 

<1 

50 

mw 

58.1 ( .5) 

54.2(1) 

<1 

4 ( . 5) 

60 

mm 

46.0(3) 

39.7(1) 

<1 

5.2(1) 

70 

mm 

28.4 ( .8) 

23.6(1) 

1-5 ( .5) 

7.3(1) 

80 

mm 

17.7(2) 

15.6(1) 

2.6 ( . 5 ) 

8 (1) 

90 

50 

12.6 ( .7 ) 

11.4 (.8) 

3 • 7 ( . 5 ) 

8.5(1) 

100 

60 

8 (1) 

7.3 (.9) 

4.4(1) 

11 (2) 


Note: Detector was located at 250 cm in free space. 
Neutrons were turned 'off* in these runs. 

Numbers in brackets are the standard deviations 
of the values they follow. T50 and T90 are the 
time intervals over which the first 50% and 90% 
(respectively) of the Dose arrives. 


In the calculations presented in Table 3.1, the inner radius 
of the LiH shell was held fixed, so that increasing the 
thickness of the LiH shell also increased the diameter. 
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The detector used in Figure 3.1 (and for most of the data 
presented in this thesis) was a point detector located 250 
cm from the center of the chamber in free space. At 250 cm, 
spherical divergence of the source pulse would lead to less 
than a 15 percent variation in fluence over a 2 m surface 
located perpendicular to the axis of the chamber. The 
uniformity of the X-ray fluence is important to the fidelity 
of the simulation, in that a space vehicle would experience 
an almost planar X-ray pulse, and many effects can be 
geometry dependant. 

Increasing the thickness of the LiH shell lengthened the 
arrival interval (T90) of the X-ray energv fluence from 3.5 
ns at a 10 cm thickness to 10 ns at a 60 cl thickness. 
Similarly, the rise-time (T50) increased from less than 1 ns 
at 10 cm to 4 ns at 60 cm. In the process, the energy 
fluence was attenuated by 90 percent of the free field 
value. 

The physics behind the relationship of the pulse length and 
rise-time to the thickness and outside diameter of the LiH 
shell can be explained as follows. When the initial pulse 
of X rays enters the inner wall of the LiH, many of them 
undergo multiple scatters before escaping through the outer 
wall of the LiH shell. The more an individual X ray 
scatters, the longer its path length and, therefore, the 
longer the time it takes to escape from LiH shell. As the 
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thickness of the shell is increased, the number of X rays 
that undergo multiple scatters increases as does the average 
number of scatters per photon; the collided pulse becomes 
longer (temporally) and takes on a longer rise-time, while 
the uncollided pulse is attenuated. Since the source pulse 
is instantaneous, the uncollided X-ray pulse has an 
instantaneous rise-time and no duration. As the ratio of 
collided X-ray fluence to uncollided X ray fluence 
increases, the pulse length and rise-time of the collided X 
rays increasingly dominate the temporal characteristics of 
the overall pulse. In the spherical geometry used in this 
work, there are two material thicknesses which affect the 
temporal qualities of the pulse: the thickness of the LiH 
shell and its outside diameter. The thickness of the shell 
determines the size of the uncollided portion of the pulse 
and, therefore, is the dominant characteristic governing the 
rise-time of the pulse. The diameter of the shell is the 
maximum thickness of LiH through which a singly-scattered X 
ray must pass, and is therefore strongly tied to the total 
pulse length. 

One measure of the duration of a pulse of X rays is the Full 
Width at Half Maximum (FWHM) of the X-ray pulse. A plot of 
the normalized X-ray flux and dose rate (rads(Si)/sec) as a 
function of time outside a 60 cm LiH shell is presented in 
Figure 3.2. The standard deviation of each point is plotted 
beneath each of the dose and flux curves. The standard 
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deviation of each point is the difference between the 
standard deviation curve and the x axis. Both the dose rate 
and the flux have a FWHM of 7 to 8 ns for a 60 cm shell. 
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Spectral Shifts 


To assess the effect of the scattering medium on the energy 
distribution of the X-ray spectrum, the resultant spectra 
from the calculations presented above were compared to the 
input X-ray spectrum. A plot comparing the input spectrum 
to the spectrum outside a 60 cm LiH shell is presented in 
Figure 3.3. LiH was chosen for the scattering medium 
because it had a very low absorption cross section for hard 
X rays. A hard X ray undergoing a collision in a LiH is 
usually scattered to a lower energy instead of being 
absorbed. On the other hand, the cross section for 
absorption of soft (low-energy) X rays is much higher than 
for hard X rays, so that a soft X ray undergoing a collision 
is almost always absorbed. When a pulse of X rays with a 
wide energy band enters the LiH shell, the high-energy X 
rays scatter down in energy and the low-energy X rays are 
absorbed. Thus, both ends of the spectrum diminish in 
magnitude, leaving a peak somewhere in the middle. The 
energy of this peak is determined by the characteristics of 
the scattering material, being higher for some than for 
others. In the LiH used for this work the peak energy 
shifted from the input spectrum peak of 37.5 keV to a final 
spectral peak of 25 keV. The final spectrum is more 
strongly peaked than the initial spectrum, an effect which 
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is more pronounced as the thickness of the LiH shell is 
increased. 


Normalized joules/eV 



1.0E-04 1.0E-05 1.0E*06 1.0E-07 1.0E-08 

Photon Energy (eV) 

10 keV Blackbody Source 

Figure 3.3: Peaking of the X-Ray Spectrum by 7 LiH 


Neutron Effects 




To assess the effects of the neutron component of the DT ICF 
pulse, a calculation was performed with the X rays turned 
'off' and the neutrons turned •on'. Neutrons interact with 
the materials in the chamber to produce neutron-capture 
gamma rays which could interfere with temporal and spectral 
fidelity of a hard X-ray simulation test. The neutrons 
themselves can be a problem if their energies and fluences 
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are sufficiently high. The gamma-ray spectrum and arrival 
times outside a 40 cm thick LiH shell in an LMF chamber are 
presented in Figure 3.4. The gamma-ray pulse consists of a 
strong peak around 2 MeV with a long, low-energy tail. The 
total gamma fluence at 250 cm was 36.6(+/- 1) cal/cm 2 and 
produced a dose of 480 (+/- 10) krad(Si). The gamma pulse 
began at 19 ns, but ninety-five percent of the pulse arrived 
after 100 ns. The total neutron fluence was 2(+/- 0.5)*10 15 
MeV-equivalent neutrons/cm 2 . An examination of the 
secondary generation data produced by MORSE revealed that 
the majority of the energy of the neutron-generated gamma 
rays was formed by the capture of epithermal neutrons by the 
lithium hydride. The large number of scatters and the 
length of time necessary for the neutrons to slow down to 
energies where the capture cross section becomes appreciable 
made the process very time-consuming. 







Figure 


Normalized Joules/ev 



Photon Eneroy (ev) 


Fraction of Total Fluence 



Time (sec) 

Note: 40..80 cm LiH Shell. Neutron Source Only 

3.4: Neutron-generated Photons Outside a 7 LiH Shell 
in an LMF Chamber 
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Combined Calculations (Neutrons + X Rays 


To examine the combined effects of the source X rays and the 
neutron-generated gamma rays, separate calculations were 
performed with the LiH shell (40 cm thick) in free space and 
in the LMF chamber. A comparison of the temporal response 
at 250 cm is presented in Figure 3.5. Twenty percent of the 
total photon energy pulse was from gamma rays generated in 
the aluminum wall, 50% was from gamma rays generated in the 
LiH and only 30% was the original X-ray pulse. The temporal 
quality of the X-ray component, however, was not actually 
changed. Ninety percent of the total photon pulse which 
arrived by 100 ns arrived within the first of 10 ns of the 
total pulse? ninety-seven percent of the gamma-ray pulse 


arrived later than 100 ns. 







Fraction of Total Fluence 



1.0E-09 1.0E-08 1.0E-07 1.0E-06 1.0E-05 1.0E-04 1.0E-03 1.0E-02 1.0E-01 

Time (sec) 

40..80 LiH Shell. 

10 keV X-Ray Source with Neutrons. 

Figure 3.5: Arrival Times for X Rays and Gamma Ray'; Outside 
a 7 LiH Shell in an LMF Chamber with Neutrons ON 


Response at a Graphite Cone 

Since most of the energy of the gamma ray pulse is carried 
by high-energy photons, it might be argued that they are 
not significant due to their small interaction probability. 
Indeed, the unfiltered (surface) dose outside a 40 cm thick 
LiH shell was 15.6 Mrad(Si) from X rays compared to 0.48 
Mrad(Si) from gamma rays. The filtered dose, however, is 
quite a different matter. After passing through the 
graphite surface of an RV and the aluminum casing of an 
electrical circuit, the high-energy component of the 
spectrum takes on a much greater significance due to the 
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attenuation of the lower-energy X rays. Figure 3.6 presents 
a comparison of the dose arrival time at a detector just 
outside the graphite RV shape in the LMF chamber to that of 
a detector located at the center of the aluminum/silicon 
shell inside the RV. Inside the graphite, the late-time 
photons comprise over 90 -cent of the total fluence while 
comprising only 70 percent at the surface. About 30 percent 
of the surface dose (rads(carbon)) arrives after 100 ns, but 
still over 90 percent of the dose (rads(C)) at depth is 
generated by late-time gamma rays. Only a few percent of 
the silicon surface dose is late-time, but at depth the 
silicon dose still exhibits the strong stepped temporal 
characteristic with 70 percent of the total dose arriving 
after 100 ns. 








Figure 3.6: Dose 














Reduction of Neutrons and Keutron->,apture Gamma Rays 

A number of methods were examined to reduce the energy 
contribution from the gamma rays and to shorten the interval 
over which they were generated. Most of the gamma rays from 
the LiH were created by the capture of epithermal neutrons 
by the lithium. The isotope of lithium used was 7 Li, which 
captures high-energy (MeV range) neutrons and low-energy 
(epithermal) neutrons, but does not capture intermediate 
energy neutrons. The addition of 10 B (which has a high 
neutron capture cross section over a broad range of 
energies) at 10% normal atom density to the LiH reduced the 
time interval over which the gamma rays were generated in 
the scattering shell to less than 100 ns. These results are 
illustrated in Figure 3.7. At 250 cm, the gamma ray fluence 
was reduced to 7.1(+/-0.6) cal/cm 2 and the gamma ray dose to 
160 (+/- 10) krads(Si). The neutron fluence was 9(+/- 
1)*10 14 Mev-equivalent neutrons/cm 2 . The 2 MeV gamma ray 
peak in the previous gamma-ray spectrum was eliminated but 
was replaced by a lower 0.478 MeV 10 B gamma-ray peak; the 
resulting photon spectrum had an endpoint energy of around 1 
Mev. The addition of boron did not affect the gamma rays 
which were generated in the walls, nor did it significantly 
reduce the MeV-equivalent neutron fluence. 
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To reduce the fluence of neutrons scattered from the wall 
and the photon fluence from gamma rays generated in the 
wall, two shielding configurations — one rectangular and 
one conical — were used. Both shields consisted of a 1 m 
external blanket of LiH seeded with 10 percent boron to stop 
the neutrons, and a 5 cm internal lead liner to stop the 
neutron-capture gamma rays. Both shields began at 150 cm 
and extended to the wall, surrounding the detector locations 
and opening towards the borated LiH shell. Runs were made 
with the 40..80 (inner radius..outer radius), 60..80, and 
20..40 cm scattering shells of borated lithium hydride. The 
temporal and spectral photon responses for both shields were 
comparable; the photon response for the 40..80 shell is 
presented in Figure 3.8. At locations deep inside the 
shield, the photon energy fluence after 100 ns (which 
includes the wall gamma rays plus late gamma rays from the 
boron) was reduced from 30% to 10% of the total photon 
energy fluence. The responses from the thinner shells were 
roughly equivalent to each other. However, they differed 
from that of the 40..80 shell in that their temporal 
response rose rapidly to a plateau at 60 percent of the 
total response by about 20 ns, then rose gradually in 
several steps and slopes to better than 98 percent of the 
total response at one microsecond. The total pulse lengths 
were still longer than desired and the MeV-equivalent- 
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neutron fluence was not significantly affected by the 
presence of the shields. 


Normalized Fluence 



Time (sec) 

40..80 LiH. X-Ray ♦ Neutron Source. 

Figure 3.8: Shielded Photon Fluence Outside a 7 LiHB Shell in 

an IMF Chamber 

Since the neutron- 7 Li reaction was the major contributor of 
gamma rays, a number of calculations were made using pure 
6 LiH as the scattering medium. These calculations were made 
with the RV turned "off" and with a cylindrical shield about 
the detector locations. The shield extended from 150 cm to 
the wall, had an inner radius of 100 cm, and had a 20 cm 
thick outer layer natural tungsten and a 20 cm thick inner 
layer of LiH. The previous chamber wall (1 m thick 
aluminum) was replaced with a more realistic wall composed 
of an 8 cm aluminum inner wall, aim thick blanket of 









borated water (0.01% 10 B), and an 8 cm aluminum outer wall. 
The fluences and doses from these calculations are presented 
in Table 3.2; the photon spectrum is presented in Figure 3.9 
and the arrival times for various thicknesses of LiH 
scattering shells are presented graphically in Figure 3.10. 
In these calculations the diameter of the LiH shell was held 
constant and the inner radius varied. As a result, the 90% 
arrival times were roughly equivalent for all thicknesses of 
LiH but the rise-time increased with the shell thickness. 


Table 3.2: Variation in Photon Fluence, Dose and Arrival 
Interval with 6 LiH Thickness 


Inner 

radius 

(cm) 

LiH 

Shell 

(cm) 

Photon 

Fluence 

(cal/cm 2 ) 

Photon 

Dose 

(MRads(Si)) 

Neutron* 
Fluence 
(10 14 /cm 2 ) 

80 

00 

80(1) 

72 (2) 



10 

51(3) 

47(2) 

300(100) 

60 

20 

36(2) 

32 (2) 

50(18) 

50 

30 

25(1) 

23 (2) 

30(9) 

40 

40 

18(1) 

18(1) 

6(2) 

30 

50 

14(2) 

13(1) 

3(0.9) 

10 

60 

10.8(1) 

9.3(.7) 

0.8(0.2) 


*Read as 10 14 Mev-equivalent neutrons/cm 2 . 

Note: Detector was located at 250 cm in free space. 

Neutrons were turned 'on' in these runs. Numbers 
in brackets are the standard deviations of the 
values they follow. 


The minimum neutron fluence in these runs was 8 (+/- 2)*10 ljl 
neutrons/ cm 2 , much closer to the desired levels, and the 
contribution of neutron capture gamma rays was negligible. 
The photon spectrum was essentially identical to the X-ray 
spectrum, with no evidence of a high-energy tail. At a 
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thickness of 60 cm there remained sufficient fluences of 
photons to be of interest to VLE work. Further, the use of 
a more realistic chamber wall seems to have eliminated the 
large neutron-capture gamma-ray contribution from the 
aluminum. 


Normalized joules/eV 



1.0E-04 I.C'E-05 1.0E-06 1.0E-07 1.0E-08 

Pnoton Energy (eV) 

Figure 3.9: Photon Spectrum (X + Gamma Rays) Outside a 60 cm 

6 LiH Shell (Neutrons ON) 
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Fraction ol Total Fluence 



Fraction of Total Photon Dose (rads(Si)) 



10 keV X-Ray Source with Neutrons. 

Figure 3.10: Temporal Variations of X-ray Dose and Fluence 
with 6 LiH Thickness (Neutrons ON) 
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A broad variety of shadow shields were used in an attempt to 
reduce the neutron fluences while retaining sufficient X-ray 
intensity to be of interest to VLE work. Without a shadow 
shield, the smallest one-MeV-equivalent neutron fluence 
calculated in this work was 8*10 13 neutrons/cm 2 . Although 
the results were not conclusive, calculations using a 
tungsten shadow shield inside the LiH shell indicate that it 
may be possible to reduce the neutron fluence by 50% while 
reducing the X-ray intensity by the same percentage. 
Reduction below 50% through the use of shadow shields alone 
seems unlikely, however, since the presence of the shield 
seems to increase the scattered contribution from the region 
of the LiH shell not shadowed by the shield. 
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IV. Summary and Recommendations 


summary 


Using a scattering medium of 6 LiH, it is possible to 
temporally stretch the initial X-ray pulse so that 90 
percent of the X-ray energy incident on the detector arrives 
over a period of 15 ns. In doing so, 87 percent of the 
initial X-ray energy (with photon energy greater than 
10 keV) is absorbed in the LiH, and the spectrum of the X- 
ray pulse is significantly softened. The remaining X-ray 
pulse retains sufficient intensity and spectral quality to 
be of interest to the VLE community. The total X-ray dose 
outside a 60 cm shell of 6 LiH was 9.3 Mrads(Si) with 50 
percent arriving over an interval of 6 ns. This corresponds 
to a dose rate 8*10 14 rads(Si)/sec. Messenger and 
Ash(2:412,413) cite hardness levels for a number of devices 
that are in the range of doses found in these calculations; 
the dose rates encountered in this work, however, exceeded 
the highest dose rates given in any of the charts presented 
in their text (10 12 rads(Si)/sec) indicating that additional 
temporal stretching and dose attenuation may be 
necessary(2:421-428). 
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The large fraction of neutrons produced by the DT reaction, 
however, pose a significant engineering problem to the use 
of the IMF as an X-ray weapons effects simulator. Neutrons 
scatter off the LiH, losing energy until they reach an 
energy level at which the neutron capture cross section is 
appreciable. When the neutrons are captured, some of the 
capturing nuclei emit high-energy gamma rays, which 
themselves scatter and contribute to the total photon 
fluence on the target. For a LiH shell containing only 7 Li, 
this process takes an extended period of time, when compared 
to the time scale of initial pulse, and could be a real 
problem if the diagnostics cannot discriminate against the 
late-time photons generated. Adding 10 B to the LiH to 
absorb the neutrons before they interact with the 7 Li does 
not completely solve the problem, but substitution of the 
isotope 6 Li for 7 Li reduces the intensity of the pulse of 
gamma rays generated in the LiH to levels which are 
insignificant when compared to the X-ray intensity. This is 
not a particularly attractive solution, however, due to the 
high cost of 6 Li. Further work needs to be done to 
determine the isotopic composition of Li which would 
maintain a sufficiently low level of 7 Li neutron-capture 
gamma rays. 
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Elimination of significant intensities of neutron-capture 
gamma rays does not eliminate the problems the neutrons pose 
to hard X-ray weapons effects simulation. The fluence of 
the neutrons themselves is high enough to present severe 
problems to tests involving semiconductor devices. In the 
calculations performed for this work, about 87 percent of 
the total source energy was in the form of neutrons. The 
minimum neutron fluence obtained in this work was on the 
order of 10 14 n/cm 2 . According to Messenger and 
Ash (2:200,401), neutron effects on semiconductors begin at 
around 10 12 n/cm 2 with many devices hardened to the 10 13 
n/cm 2 range. Calculations using shadow shields, though 
inconclusive at this point, indicate that it may be possible 
to reduce the neutron fluence by 50% and still retain 
sufficient X-ray intensity. If a suitable shielding 
configuration is discovered the neutron fluences may not be 
a problem. If real-time diagnostics are used in the 
simulations, it may be possible to diagnostically separate 
the neutron effects from the X-ray effects. 14 MeV neutrons 
travel much slower than photons; the neutrons arrive at a 
point 325 cm from the source 60 ns after the photons. If 
the diagnostics can take advantage of this interval, 
meaningful X-ray information could be gained even if the 
neutrons interfered with the effects being studied. 
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Recommendations for Further Work 


This work examined a large number of chamber configurations 
using LiH as a scattering medium, and attempted to minimize 
the fluence of gamma-rays generated by the neutrons from the 
DT reaction. This effort included adding shields and 
reflectors, many of which are not discussed here, but all of 
which had severely limited usefulness. 10 B was added to the 
LiH to eliminate the problem of gamma rays generated in the 
LiH, and while this did reduce the gamma intensity and pulse 
length, it was not as effective as could be hoped. The 
substitution of 100% 6 Li for the lithium in the LiH was 
found to reduce the gamma ray intensity to sufficiently low 
levels for VLE work. Further work is needed to determine 
the exact isotopic composition of the lithium required to 
keep the gamma ray levels sufficiently low. Other 
scattering media, such as carbon and beryllium, as well as 
other neutron absorbers such as cadmium, hafnium and indium 
should also be investigated. Boron was an interesting 
absorber choice because its use did not introduce an 
activation problem. Other choices will need to be carefully 
evaluated in this regard. 

The use of shadow shields offers great promise in further 
reducing the fluence of energetic neutrons. The most 
promising shield was a conical shield composed of a layer of 
tungsten (to slow down the fast neutrons) followed by a 
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dense borated hydrogenous medium (to further slow and absorb 
the neutrons), followed by another tungsten layer (to stop 
the gamma rays). 
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Conclusion 


While the neutrons produced in the DT reaction pose 
significant problems to the use of the LMF to simulate the 
effects of hard X rays from nuclear weapons, there is 
evidence that these problems are not insurmountable. 
Although further work in this area is needed, adequate 
shielding should be able to reduce the neutron fluence to 
acceptable levels while stretching the temporal response of 
the X-ray pulse to a length, intensity and spectral quality 
to be of interest to the VLE community. 
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Appendix A 

Pser-vritten Routines 


Overview 

The MORSE computer program is supplied as a series of 
Fortran subroutines that perform the functions common to all 
problems the program is designed to handle. The user 
supplies the problem-specific subroutines, and compiles and 
links all the subroutines together to form the executable 
program. This appendix presents a brief discussion of the 
function of each subroutine used in this work, followed by 
the source code for all the subroutines. The subroutines 
were written in Fortran 77 and were compiled and run on a 
VAX under the VMS operating system. 
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Description of the Subroutines 

Main . This is the main controlling module of the 
program. It allocates all of the common memory blocks, 
reads in the input and output file names, and terminates the 
program. 

Bankr . Bankr is, as its name implies, the banker of 
the MORSE system. Just like an actual banker, it is 
consulted any time a significant event occurs in the life of 
a particle, such as a birth, collision or boundary crossing. 
It is Bankr that calls the flux estimator routines Relcol, 
Sdata and Sgam. 

Direc . Direc is a required user-written function. It 
must be in the code, but doesn't really do anything in this 
application. 

Gtmed . The purpose of this function is to assign a 
different media to a particular zone other than that 
assigned in the data file. In this application, the routine 
simply sets the media for each zone to that assigned in the 
data file. 

Relcol, Sdata and Sgam . Relcol is a routine used to 
estimate the flux at a point detector. At each collision, 
Bankr calls Relcol which determines (through subroutine 
calls) the probability per unit steradian of a particle 
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being emitted in the direction of the detector, the number 
of mean free paths to the detector, and the distance to the 
detector. 
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The contribution of the particle to the total fluence at the 
detector is then calculated as (5:4.6-39): 

Con=Wate*Pro(mu,E)*exp(-mfp)/Dist 2 , (1) 

where: 

Con is the contribution of the particle to the 

fluence at the detector, 

Wate is the weight of the particle leaving the 

collision site, 

Pro(mu,E) is the probability per unit steradian that 

the particle scatters to outgoing energy 
group E with direction cosine mu, 

exp() is the exponential operator, 

mfp is the number of mean free paths between the 

collision site and the detector, and 

Dist 2 is the square of the distance to the 

detector. 

This contribution is then multiplied by the average energy 
of the particle energy bin (taken as the simple average of 
the upper and lower energy of the bin) to produce an energy 
fluence in ev/cm 2 , and passed to the bookkeeping routines 
along with the age of the particle. 



It should be noted, however, that Eq.(1) is only good for 
contributing events that occur at distances greater than one 
centimeter from a detector. Nearer than this the formula 
begins to blow up, approaching infinity as the distance 
approaches zero. For this reason, detectors should always 
be located one centimeter from the nearest surface on which 
a source event or collision could occur; the code contains a 
provision to skip this calculation and print a warning 
should the distance to the detector be less than one 
centimeter. 

Sdata and Sgam perform the same function as Relcol (and are 
actually modifications of Relcol), but at different points 
in the particle lifetime. Sdata sends information to the 
bookkeeping routines when the particle is generated at the 
beginning of the batch. These data form the uncollided 
portion of the total fluence, and the probability per 
steradian of a particle having direction cosine mu is taken 
as 1/(4 pi). Sgam sends information at the birth of a 
secondary particle, and is identical to Sdata except that 
the fluence information is stored as collided fluence. 

There are a number of TYPE statements included in these 
routines. These are to print diagnostic information in case 
the program encounters a geometry problem. 
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Source . Source is the subroutine that generates the 
source particles. It determines the energy group of the 
source particles from the input energy spectrum and 
determines their direction cosines randomly from a discrete 
isotropic angular distribution. Source also sets the 
initial weight of each particle to the weight specified by 
the user in the data input file. 








Fortran Source Code 


c * 
c * 
c * 
c • 

c 

c«* 

c* 

c* 

c* 

C**’ 

c 

c * 
c * 
c 

c 

c * 
c * 

c * 
c * 


c * 
c * 


c * 
c * 


c * 
c * 


c * 
c * 

c * 


100 


200 


300 


* THIS IS THE KA!H ROUTINE •••*•••*•*•**•* ‘ * 

* Based on a routine by 0. Belter Aug 8fi. 3 blanket regions, 

* 4 detectors 

* Modified by Don Nichols Sept SB for RV in LNF chanter 


This version determines uncollided fluenee using estimators 
to calculate the contribution of each source and collision 
event to the flux at a point detector. 

******«****»******«r*ft***+******##*****'»****»***************l 


* THE FOLLOWING CARD DETERMINES THE SIZE ALLOWED FOR BLANK COMMON 

* The value of NLFT below should be set to one less than this size 


COMMON NCC600001) 


* NOTE - THE ORDER OF COMMONS IN THIS ROUTINE IS IMPORTANT AND MUST 

* CORRESPOND TO THE ORDER USED IN DUMP ROUTINES SUCH AS KELP, XSCHLP 

* 

* LABELLED COMMONS FOR WALK ROUTINES ••******•**•*• 
COMMON /APOLLO/ AGSTRT,DDF.DEADWT(26),I TOUT,IT1N 

COMMON /FISBNK/ MFISTP 
COMMON /NUTRON/ NAME 

* 

* LABELLED COMMONS FOR CROSS-SECTION ROUTINES *«•*•*•*•*• 
COMMON /LOCSIG/ ISCCOG 

COMMON /MEANS/ NM 
COMMON /MOMENT/ NMOM 
COMMON /QAL/ O 
COMMON /RESULT/ POINT 

* 

* LABELLED COMMONS FOR GEOMETRY INTERFACE ROUTINES • *«*••*• 
COMMON /GEOMC/ XTWO 

COMMON /NORMAL/ UNORM 

* 

* LABELLED COMMONS FOR USER ROUTINES • *••»••*••••*•* 
COMMON /PDET/ NO 

COMMON /USER/ AGST 

* 

* COMMON /DUMMY/ WILL NOT BE FOUND ELSEWHERE IN THE PROGRAM * * * * 
COMMON /DUMMY/ DUM 


CHARACTER*40, NAH1 
CHARACTER*4C, NAM2 
TYPE 

********* MORSE Code, LMF-ICF Problem *•••»**•*> 
=««="> WARNING !!! <===*====' 

ABORT if mixed x-secs are not assigned to FOR010 1 


TYPE 

TYPE 

TYPE 

TYPE 

TYPE 


'ENTER NAME OF INPUT FILE' 
ACCEPT 100, NAM1 
FORMAT(A40) 

TYPE VENTER NAME OF OUTPUT FILE' 
ACCEPT 200, NAM2 
FORMAT (A40) 

OPEN(UNIT=1,NAME=NAM1,TYPE='OLD') 
OPEN(UN IT=2,NAME=NAM2,TYPE 1 'NEW') 
ITOUT * 2 
1TIN « 1 
NLFT=600000 
CALL MORSE (NLFT) 

TYPE 300, NAM2 

FORMAT(X,'OUTPUT FILE IS ',A40) 
STOP 
END 
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SUBROUTINE BANKR(NBNKID) 

C DO NOT CALL EUCLID FROM BANKR(7) 

COmOH /APOLLO/ AGSTRT,DDF,DEADWT(5),ETA,ETATH,ETAUSD,UINP,V1NP, 

1 W1NP,WTSTRT,XSTRT,YSTRT,ZSTRT,TCUT,XTRA(10), 

2 10,II,MED IA,1ADJM,ISBI AS,I SOUR,ITERS,IT IME,ITSTR,LOCVTS,LOCFWL, 

3 LOCEPR.LOCNSC,LOCFSN,MAXGP,MAXTIM.MEDALB,MGPREC,MXREC,NALB, 

4 NDEAD(5),N£WNM,NGEOM,NGPOT1,NGPOT2,NGPOT3,NGPOTG,NGPCTN,NITS, 

5 NKCALC,NX ILL,NLAST,NMEN,NMGP,NMOST,NMTG,NOLEAK,NORMF,NPAST, 

6 NPSCL(13),NOUIT,NSIGL,NSOUR,NSPLT,NSTRT,NXTRA(10) 

COMMON /NUTRON/ NAME.NAMEX,IG,IGO,NMED,MEDOLD,NREG,U,V.U.UOLD,VOLD 
1 ,WOLD,X,Y,2,XOLD,YOLD,ZOLD,UATE,OLDWT,WTBC,BL2NT,BL ZON,AGE,OLDAGE 
NBNK * NBNKID 
IF (NBNK) 100,100,140 

100 NBNK * NBNK ♦ 5 

GO TO (104,103.102,101),NBNK 

101 CALL STRUN 

C CALL HELP(4HSTRU,1,1,1,1) 

RETURN 

102 NBAT = NITS - ITERS 
NSAVE * NHEN 

CALL STBTCH(NBAT) 

C NBAT IS THE BATCH NO. LESS ONE 

c This lets us know where we are in the run 

TYPE 3333, NBAT+1 

3333 FORMAT (IX.'Running batch Number',2X,14) 

RETURN 

103 CALL NBATCH(NSAVE) 

C NSAVE IS THE NO. OF PARTICLES STARTED IN THE LAST BATCH 
RETURN 

104 CALL NRUN(NITS.NOUIT) 

C NITS IS THE NO. OF BATCHES COMPLETED IN THE RUN JUST COMPLETED 
C NOUIT .GT. 1 IF MORE RUNS REMAIN 

C .EO. 1 IF THE LAST SCHEDULED RUN HAS BEEN COMPLETED 

C IS THE NEGATIVE OF THE NO. OF COMPLETE RUNS, WHEN AN 

C EXECUTION TIME KILL OCCURS 

RETURN 


140 GO TO (1,2,3.4,5,6,7,8,9,10,11,12,13),NBNK 


c 

NBNKID 

COLL TYPE 

BANKR CALL 

NBNKID 

COLL TYPE 

BANKR CALL 

c 

1 

SOURCE 

YES 

(HSOUR) 

2 

SPLIT 

NO 

(TESTW) 

c 

3 

FISSION 

YES 

(FPROB) 

4 

GAMGEN 

YES 

(GSTORE 

c 

5 

REAL COLL 

YES 

(MORSE) 

6 

ALBEDO 

YES 

(MORSE) 

c 

7 

BORYX 

YES 

(NXTCOl) 

8 

ESCAPE 

YES 

(NXTCOL 

c 

9 

E-CUT 

NO 

(MORSE) 

10 

TIME KILL 

NO 

(MORSE) 

c 

11 

R R KILL 

NO 

(TESTW) 

12 

R R SURV 

NO 

(TESTW) 

c 

13 

GAMLOST 

NO 

(GSTORE) 






1 CALL 

SDATA 








RETURN 








2 RETURN 








3 RETURN 








4 CALL 

SCAM 








TYPE 

*,'another 

gamma 

in batch 1 

,-.baf1 





RETURN 

5 CALL RELCOL 
nuncol =nuncol*1 

if (fx/ncol .GT.50Q) then 
ncncol^O 

TYPE *,'another 500 collisions in batch ',nbat + 1 
endi f 
RETURN 

6 RETURN 

7 continue 

c TYPE ‘.'another boundary crossing in batch ‘,nbat+1 
RETURN 

8 RETURN 

9 RETURN 

10 RETURN 

11 TYPE *,'another one bites the dust in batch ‘,rfcat+1 
RETURN 

12 RETURN 
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13 RETURN 
END 
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FUNCTION DIREC(F) 

direc*1. 

return 


end 

C * • 

C * * 

C * * Subroutine GTMED Bets the media for each zone 
C this version ';<)ch zone has it's own cross section assiflned 
SUBROUTINE GTH'y(MOGEOH,HDXSEC) 

MDXSEC=HDGEOM 

RETURN 

END 
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SUBROUTINE RELCOL 
C 

C THIS VERSION IS FOR POINT DETECTORS LOCATED AT (XD,YD,ZD) 

C 

COMMON /USER/ AGSTRT,UTSTRT,XSTRT,YSTRT,ZSTRT,DFF,EBOTN,EBOTG, 

1 TCUT,10,11,IADJM.NGPQT1,NGPQT2,NGPOT3,NGPQTG.NGPQTN,NITS,NLAST, 

2 NL£FT,NMGP,NMTG,NSTRT 

COMMON /POET/ ND,NNE,NE,NT,NA,NRESP,NEX,NEXND,NEND,NDNR,NTNR,NTNE, 

1 NANE.NTNDNR,NTNEND,NANEND,LOCRSP,LOCXD,LOCIB,LOCCO,LOCT,LOCUO, 

2 LOCSD,LOCQE,LOCOT,LOCQTE,LOCQAE,LMAX,EF1RST,EGTOP 

COMMON /NUTRON/ NAME,NAMEX,IG,IGO.NMED,MEDOLD,NREG.U,V.W.UOLD,VOLD 
1 ,UOLD,X,T,Z ( XOLD, TOLD, ZOLD,MATE,OLDWT ,UTBC,BLZNT,BLZON, AGE, OLDAGE 
COMMON BL(1) 

DIMENSION NL(1) 

EQUIVALENCE (BL(1),NL(1)) 

DATA NEST /I/, FNEST /I./ 

C NEST ♦ FNEST ARE THE NO. OF ESTIMATES TO BE MADE TO EACH DETECTOR 
C * * * ISTAT MUST BE EOUAL TO 1. 

C * * * NEX MUST BE AT LEAST 1 
C * * * NEXND MUST BE AT LEAST 1 
C ... 

C THIS ROUTINE REDUCES THE NUMBER OF THERMAL NEUTRONS TRACKED 

C THRU THE GEOMETRY BY 90X AND THEN COMPENSATES BY 

C MULTIPLYING THE WEIGHT OF THE REMAINING 10X BY 10. 

c ................................................................ 

CHECKS. 

IF (IGO.NE.28) GO TO 100 
Xl = FLTRNF(O) 

IF (XI.GE.0.1) RETURN 
CHECK=10. 

100 CONTINUE 

DO 30 1=1,ND 
IA=LOCXD+I 
XE = BL(IA) 

YE = BL(IA.NO) 

ZE * BL(IA*2*ND) 

A = XE • X 

B = YE - Y 

C = ZE - Z 

SD2=A*A+B*B+C*C 
DS=SQRT (SD2) 

C * » • COS DEPENDS ON THE ANGLE OF INTEREST 
COS=C/DS 

THETA = (A*UOLD + B*VOLD ♦ C*WOLD)/DS 
I GOLD = ICO 
I GO = NGP0T3 

IF (IGO.LE.NGPQT1) IGO=NGPQT1 

IA = LOCRSP ♦ NRESP*NMTG ♦ 1 

CALL PTHETA(NMED, I GOLD,! GO,THETA,BKIA), NMTG) 

NES = 0 
PSUM = 0. 

IA = IA - 1 
DO 5 IL = I GOLD,I GO 
5 PSUM = PSUM ♦ ABS (SLCIA+IL)) 

10 R = FLTRNF(0) * PSUM 
DO 15 IL = I GOLD,I GO 
IF (R-ABS (BL(IA+IL))) 20,20,15 
15 R = R - ABS CBL<I A*IL>> 

IL = I GO 
20 MARK=1 

AGED * AGE ♦ DS/BL<NMTG*IL) 

MEDim=NMED 

CALL EUCLID(MARK,X,Y,Z,XE,YE,ZE,DS,IL,ARG,0,MEDIUM,BLZNT,NREG) 
if (ARG.NE.-160) GO TO 21 
TYPE *,*g1 IT HAPENED AGAIN' 

TYPE 300,1 

300 FORHATJX, 'DETECTOR ',6X,K) 

TYPE 301, NMED 

301 FORMAT(X,'-MEDIUM ',6X,I4) 
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TYPE 302, MEDOLD 

302 FORMAT(X, 1 -OLDNED1UM '.PX.K) 

TYPE 303, NREG 

303 FORMAT(X,'-REGION *,6X,I4) 

TYPE 304, IBIZN 

304 FORMATCX,'-ZONE ',6X,I4) 

TYPE 30S.IBLZO 

305 FORMAT(X,•-OLDZONE ',6X,14) 

TYPE *, 'RELCOL' 

21 IF (ARG.LT.-64.) GO TO 25 

C ........ BEUARE TH!S VERSION WILL NOT WORK IF ENERGY BIASING IS USED 

IF (S02.LT.1.0) then 

TYPE *, 'Thre* out a particle...WARNING.* 

GO TO 25 
endif 

CON « CHECK*WATE*EXP (ARG)*SIGN (PSUM,BL(1A*1L>>/S02/FNEST 
IF (CON.IT.1.£-38) GO TO 25 
IF ((IL.NE.NGPQT2).AND.(IL.NE.NHTG)) THEN 
AVENG*(BL(I L)*BL(1L*1))/2.0 
ELSEIF (IL.EQ.NGPOT2) THEN 
AVENG=(8L(IL;*E80TN)/2.0 
ELSE 

AVENG*(BL(1L)*EBOTG)/2.0 

ENDIF 

CON=CON*AVENG 

CALL FLUXST (I,IL,CON,AGED,COS,0) 

25 NES * NES ♦ 1 

INN=LOCXD+6*ND+I 
NL(INN)=NL(INN)+1 
IF (NES-NEST) 10,30,30 
30 CONTINUE 
RETURN 
END 
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SUBROUTINE SDATA 
C 

C THIS VERSION IS FOR POINT DETECTORS LOCATED AT (XD,YD,ZD) 

C 

COMMON /USER/ AGITRT,VTSTRT,XSTRT,TSTRT,ZSTRT,DFF,EBOTN,EBOTG, 

1 TOUT,10,11,1 ADJM,NGPGT1,NGPOT2,NGPOT3,NGPQTG,NGPQTN,NITS,NLAST, 

2 NLEFT,NMGP,NHTG,NSTRT 

C0H40N /POET/ ND,NNE,NE,NT,NA,NRESP,NEX,NEXND,NEND,NDNR,NTNR,NTNE, 

1 NANE.NTNDNR,NTNEND.NANEND,LOCRSP,LOCXD,LOCIB,LOCCO,LOCT,LOCUO, 

2 LOCSO,LOCOE,LOCQT,LOCOTE,LOCQAE,LKAX,EF1RST,EGTOP 

COMMON /NUTRON/ NAME,NAMEX,IC ( I GO,NMED.MEDOLD,NREG.U,V,W.UOLD,VOLD 
1 ,WOLD,X,T,2,XOLD,TOLD,ZOLD,UATE,OLDWT,WTBC,BLZNT,BL20N,AGE,OLDAGE 
COMMON EN(1) 

FNEST*1.0 
DO 25 1*1,ND 
1A*LOCXD+! 

XE * EN(I A) 

YE * EN(IA*ND) 

ZE * EN(IA*2*ND) 

A * XE - X 
B = YE - Y 
C = ZE - Z 
SD2*A*A+B*B + C*C 
DS*SQRT (SD2) 

C * * * COS DEPENDS ON THE ANGLE OF INTEREST 
COS=C/DS 

AGED=DS/EN(NMTG*IG)+AGE 

20 MARK=1 
MEDIUM*NMED 

CALL EUCLID(MARK,X,Y,Z,XE,YE,ZE,DS,IG ( ARG,0,HEDIUM,BLZNT,NREG) 
if (ARC.NE.-160) GO TO 21 
TYPE *,'gl IT HAPENEO AGAIN' 

TYPE 300,1 

300 FORMATCX, 'DETECTOR \6X,I4) 

TYPE 301, NMED 

301 FORMATCX,'-MEDIUM ',6X,I4) 

TYPE 302, MEDOLD 

302 FORHAT(X,'-OLDMEDIUM ',6X,I4) 

TYPE 303, NREG 

303 FORMATCX,'-REGION ',6X,I4) 

TYPE 304, IBLZN 

304 FORMATCX,'-ZONE ',6X,14) 

TYPE 305,IBLZO 

305 FORMATCX,'-OLDZONE ',6X,14) 

TYPE *, 'SDATA' 

21 IF CARG.LT.-64.) GO TO 25 

C’*’*****BE U ARE THIS VERSION WILL NOT WORX IGNORES COLLISIONS CLOSE TO 
DETECTOR 

IF (SD2.LT.1.0) THEN 

TYPE *, 'Threw out a particle...WARNING.' 

GO TO 25 
END IF 

CON = WATE*EXP (ARG)/12.56637/SD2/FNEST 
IF (CON.LT.1.E-38) GO TO 25 
IF (CIG.NE.NGPQT2).AND.(IG.NE.NMTG)) THEN 
AVENG=(EN(IG)*EN(IG*1))/2.0 
ELSEIF (IG.EQ.NGPOT2) THEN 
AVENG=(EN( IG)*EBOTN)/2.0 
ELSE 

AVENG*(EN(IG)*EBOTG)/2.0 

END IF 

CON*CON*AVENG 

C * • SWITCH=1 -- STORE IN ALL RELEVANT ARRAYS 

CALL FLUXST (I,IG.CON,AGED,COS,1) 

25 CONTINUE 
RETURN 
END 
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SUBROUTINE SGAM 
C 

C THIS VERSION IS FOR POINT DETECTORS LOCATED AT CXD.YD.ZD) 

C 

CO#«ON /USER/ AGSTRT,WTSTRT,XSTRT,YSTRT,ZSTRT,DFF,EBOTN,EBOTG, 

1 TCUT,IO,11,1ADJN,NGPOT1,NGPGT2,NGPCT3,NGPGTG,NGPOTN,NITS,NLAST, 

2 NLEFT,NMGP,NMTG,NSTRT 

COFMON /POET/ NO ( NNE,NE,NT,NA,NRESP,NEX,NEXND,NEND,NDNR,NTNR,NTNE, 

1 NANE,NTNDNR,NTNEND ( NANENO,LOCRSP,LOCXD,LOCIB.LOCCO,LOOT,LOCUO, 

2 LOCSO,LOCOE,LOCOT,LOCOTE,LOCQAE,LHAX,E FIRST.EGTOP 

C0FR4ON /NUTRON/ NAME,NAMEX,IG,IGO,NM£D,MEDOLD,NREG.U,V,W.UOLD,VOLD 
1 ,WOLD,X,Y,Z,XOLD,YOLD,ZOLD,MATE,OLDWT,WTBC,BLZNT,BLZON,AGE,OLDAGE 
COMMON EN(1) 

c ................................................................ 

FNEST-1.0 
DO 25 1*1,HD 
IA*LOCXD+l 
XE « EM CIA) 

YE * ENOA+ND) 

ZE * EN(IA+2*ND) 

A = XE - X 
B * YE - Y 
C * ZE - Z 
SD2=A*A*B*B*C*C 
DS'SORT (SD2) 

C * * * COS DEPENDS ON THE ANGLE OF INTEREST 
COS=C/DS 

AGED=DS/EN(NMTG*IG)+AGE 
20 MARK=1 

NEDIUM=NMED 

CALL EUCLIDCMARK,X,Y,Z,XE,YE,ZE,DS,1G,ARG,0,MEDIUM,BLZNT,NREG) 
if (ARC.NE.-160) GO TO 21 
TYPE *,*0' IT HAPENEO AGAIN 1 
TYPE 300,1 

300 FORMATCX,’DETECTOR *,6X,I4) 

TYPE 301, NMED 

301 FORMATCX,'-MEDIUM *,6X,I4) 

TYPE 302, MEDC.D 

302 FORMATCX,'-OLDMEDIUM ',6X,I4) 

TYPE 303, NREG 

303 FORMAT(X,'-REGION ',6X,14) 

TYPE 304, IBLZN 

304 FORMATCX,'-ZONE \6X,!4) 

TYPE 305,IBLZO 

305 FORMATCX,'-OLDZONE ',6X,I4) 

TYPE *, 'SGAM' 

21 IF CARG.LT.-64.) GO TO 25 


C**.‘BEUARE THIS VERSION WILL NOT WORK IF COLLISION IS CLOSE TO 

c ........ DETECTOR 


IF CS02.LT.1.0) THEN 

TYPE *, 'Threw out a part icle...WARNING.' 

GO TO 25 

ENDIF 

CON = WATE*EXP (ARO/12.56637/SD2/FNEST 
IF CCON.LT.1.E-38) GO TO 25 
IF CCIG.NE.NGPOT2).AND.CIG.NE.NMTG)) THEN 
AVENG=CENCIG)*£NC1G*1))/2.0 
ELSE!F CIG.EQ.NGP0T2) THEN 
AVENG=CENCIG)«EBOTN)/2.0 
ELSE 

AVENG*CENC!G)<-EBOTG)/2.0 

ENDIF 

CON*CON*AVENG 

C * • SWITCH*0 -- STORE IN ALL RELEVANT ARRAYS EXCEPT UD 

CALL FLUXST (I,IG,CON,AGED,COS,0) 

25 CONTINUE 
RETURN 
END 
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SUBROUTINE SOURCE(IG,U,V,W,X,Y,Z,WATE,MED,AG,I SOUR,ITSTR,NGP0T3, 
1 DOF,I SB 1 AS,NUTG) 

COWKM /USER/ DUM(9), 10,11,1DUH( 12) 

COMMON VTS(1) 

c 

C IF ITSTR«0, MUST PROVIDE IG,X,Y,Z,U,V,W,WATE AND AG IF DESIRED TO BE 
C DIFFERENT FROM CARD VALUES (WHICH ARE THE VALUES INPUT TO SOURCE) 

C IF ITSTR*1, IG IS THE GRP NO. CAUSING FISSION, MUST PROVIDE NEW IG 

C THIS VERSION OF SOURCE SELECTS INITIAL GROUP FROM THE INPUT SPEC 
C 

DATA 1 CALL/1/ 

IF (ICALL) 10,10,5 
5 ICALL = 0 

WRITE (10,1000) 

1000 FORMAT (> YOU ARE USING THE DEFAULT VERSION OF SOURCE WHICH SETS W 
1ATE TO DDF AND PROVIDES AN ENERGY IG.') 

10 IF(I SOUR)15,15,60 

15 WATE=DDF 

IF (ISBIAS) 20,20,25 
20 NWT * 2*NNTG 
GO TO 30 

25 NWT * 3*NHTG 
30 R = FLTRNF(0) 

DO 35 1=1,NGPOT3 
IF (R - WTS(I+NWT)) 40,40,35 
35 CONTINUE 
40 IG=I 

IF (ISBIAS) 60,60,45 
45 IF (1-1) 60,50,55 

50 WATE = WATE*WTS(2*NMTG+1)/WTS(3*NMT G-*1) 

GO TO 60 

55 WATE = WATE*(WTS(2*NMTG*I)-WTS(2*NMTG*I-1))/(WTS(3*NMTG*I)-WTS(3*N 
1MTG+I-1)) 

c ****•**••*•***••*#*•»**•**•**•*••**•*******•***•*••****••**««*•* 

C This version of the MORSE SOURCE file is modified for a 

C problem with the point source directed ISOTROPICALLY. 

C *•*•*•*•*•******•****•****•••••****•«*******••*•*•*•*••*******•* 

60 Z = 0.0 

Y = 0.0 

X = 0.0 

99 CALL GTISO(U,V,W) 

RETURN 

END 
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Appendix B 


Ex ample of MORSE Data Input File 

Data is entered into the MORSE program at run time in the 
form of a data input file, with the format of each record 
precisely determined by the format statements in the source 
code and documented in the manuals that accompany MORSE 
(5 : ch 4). This appendix presents a sample data input file, 
that was used to calculate the temporal and spectral photon 
fluence at positions inside a shielded cavity in an LMF 
chamber, outside a LiHP shell. 
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